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1 Introduction 

One of the great challenges in theoretical elementary particle physics still is the description 
of bound states in a way fully consistent with all requirements imposed by special relativity 
and quantum mechanics, that is, within the framework of relativistic quantum field theory. 
In principle, the appropriate tool to achieve this goal is the Bethe-Salpeter formalism [1-3]: 
for bound-state constituents with fixed features, the homogeneous Bethe-Salpeter equation 
governs all the bound states. However, for various practical reasons some of which we recall 
below, in its applications to both quantum electrodynamics and quantum chromodynamics 
frequently some simplified equations, situated along a path of nonrelativistic reduction, are 
used. Its cornerstones are the Salpeter equation [4] and the reduced Salpeter equation [5,6]. 

The nature of such bound-state equations renders difficult to find, for given interactions 
of the bound-state constituents, their exact analytic solutions. However, a knowledge of, at 
least, some exact solutions facilitates to judge the significance of the outcomes of reduction 
steps or numerical solution techniques. In view of this, we remember a straightforward idea 
used in quantum physics, e.g., by Neumann and Wigner [7] for studies of the following kind. 

Within the formalism of nonrelativistic quantum mechanics, bound states are described 
by the time-independent Schrodinger equation, the eigenvalue equation of the Hamiltonian 
operator controlling the dynamics of the system under consideration for energy eigenvalues 
E and associated eigenf unctions ip. The bound states correspond to the discrete eigenvalues 
in the spectrum of this operator. For a particle of mass m experiencing interactions induced 
by some potential V, this Schrodinger equation in configuration-space representation reads 

A 1 

ip(x) = Eip(x) , A = V • V 



2m 



+ V(x) 



Typically, the interaction potential V is inferred from physical considerations or principles, 
and, for this potential, the solutions {(E, ■?/>)} of the Schrodinger equation then are derived. 
However, for the purpose of just constructing examples for exact solutions (E, vp) related to 
potentials V entering into the Schrodinger equation one may also follow the opposite route: 
on an equal footing, one may postulate, for a chosen value of E, any preferred or convenient 
shape of the solution ip and see whether one is able to find analytically the related potential 
V. This is easily done for vanishing binding energy E, i.e., for E = 0, with the general result 

v(x) = . 

Assuming the potential V to be spherically symmetric, V(x) = V(r), r = \x\, the choice [7] 

, , . , , . sin(a 3 r 3 ) 
ip(x) = ip(r) = — , a = 2m, 

for a likewise spherically symmetric stationary solution thus fixes [7] the central potential: 1 

V(r) = —= -9a 5 r 4 . 
ar 1 

We intend to adapt the above procedure, in quantum mechanics applied to Schrodinger 
problems, to the Bethe-Salpeter approach. To this end, we first identify, in Sec. 2, examples 
of Bethe-Salpeter equations with internal structure sufficiently simple to allow for the kind 
of inversion we have in mind. For these tractable bound-state equations, we then, in Sec. 3, 
relate postulated candidate solutions to confining and non-confining interaction potentials. 



1 Here, the eigenvalue E = is located in the continuous spectrum of the Hamiltonian, thus not discrete. 
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2 Bethe-Salpeter Formalism in Instantaneous Limit 

Motivated by the needs of physics, the aim of the present discussion is to construct analytic 
solutions to equations of motion describing relativistic bound states which we assume to be 
composed of some fermion, with momentum pi, and some antifermion, with momentum p 2 ; 
the total momentum of the bound state is P = P1+P2 and its mass M is given by M 2 = P 2 . 



2.1 Homogeneous Bethe— Salpeter Equation [1—3] 

In principle, within the framework of relativistic quantum field theory the adequate tool for 
the (Poincare-covariant) description of bound states is the Bethe-Salpeter formalism [1-3]. 
In this approach a bound state B(P) of momentum P and mass M = \f~P 2 is described by a 
Bethe-Salpeter amplitude defined in configuration-space representation as matrix element 
of the time-ordered product of the field operators of the bound-state constituents evaluated 
between vacuum state |0) and bound state |B(P)). In momentum-space representation, the 
Bethe-Salpeter amplitude, after splitting off the center-of-momentum motion of the bound 
state and suppressing all indices generically denoted by P), encodes the distribution of 
the relative momentum p of the two bound-state constituents; it satisfies the formally exact 
homogeneous Bethe-Salpeter equation, which involves two kinds of dynamical ingredients: 
the full propagators Si(pi) of the constituents with individual momenta Pi,i — 1,2, and the 
Bethe-Salpeter interaction kernel K(p, q, P), by construction a fully amputated four-point 
Green function of the bound-state constituents defined perturbatively by summation of the 
countable infinity of all Bethe-Salpeter-irreducible Feynman diagrams for two-particle into 
two-particle scattering. Skipping all indices, this Bethe-Salpeter equation generically reads 

Hp,P) = j^y A S 1 (p 1 ) f d 4 qK(p,q,P)<S>(q,P)S 2 (-p 2 ) . (1) 

Unfortunately, attempts to apply the Bethe-Salpeter formalism to relativistic bound-state 
problems are doomed to face various grave obstacles, such as the impossibility to determine 
the Bethe-Salpeter interaction kernel beyond the tight limits of perturbation theory, or the 
appearance of excitations in the relative time variable of the bound-state constituents, that 
is, of solutions called abnormal, difficult to interpret in the framework of quantum physics. 2 
They may be remedied in three-dimensional reductions of the Bethe-Salpeter equation; the 
most well-known of the resulting bound-state equations is the one proposed by Salpeter [4]. 



2.2 (Full) Salpeter Equation [4] 

In order to facilitate the formulation of instantaneous Bethe-Salpeter equations, let's recall 
the standard definition of one-particle energy Ei(p), one-particle Dirac Hamiltonian H^p), 
and energy projection operators Af(p) for positive and negative energy of particle i = 1,2: 

E z (p) ee sjp 2 + m 2 , H t (p) ee 7o ( 7 ■ p + m t ) , Af (p) ee ^f^g^ ■ 

The projection operators A.f(p) satisfy, in terms of the charge-conjugation Dirac matrix C, 
[Af(p)] c ee [C- 1 Af (p) cf = Af(p) , i = 1,2 . 



There are, however, very sound arguments [8] claiming that the presence of time- like excitations might 
be just an artifact of the ladder approximation to the Bethe-Salpeter kernel, incorporating the interactions 
responsible for the bound states only by an iteration of a single-particle exchange between the constituents. 
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The reduction of the Bethe-Salpeter equation to the Salpeter equation relies on exactly 
two simplifying assumptions (brief reviews of such reductions can be found in Refs. [9-11]): 

First, the instantaneous approximation assumes that in the center-of- momentum frame 
of the bound states, fixed by P = (M, 0), all interactions between bound-state constituents 
are instantaneous. In this ("static") limit, the Bethe-Salpeter interaction kernel K(p, q, P) 
depends only on the spatial components, p and q, of the relative momenta p and q involved, 
that is, takes the form K(p, q, P) = K{p, q). If, furthermore, the propagators Si(pi) of both 
constituents may be assumed to be entirely free from any nontrivial dependence on the zero 
component p of the relative momentum p, integrating w.r.t. p reduces the Bethe-Salpeter 
equation (1) to a kind of instantaneous Bethe-Salpeter equation for the Salpeter amplitude 

4>{P) = ^J d Po®(p) ■ 

For various attempts in this direction, see, for instance, Refs. [12,13] and references therein. 

Second, any bound-state constituent is assumed to propagate as free particle with some 
effective mass required to encompass appropriately all dynamical self-energy effects. 3 In 
other words, any fermion propagator Si(p) is approximated by its corresponding free form 4 

Sifi{p,mi) = - — = i— T —- , eiO, z = l,2. 

p — rrii + i e p z — + i e 

Imposing both restrictions above on the Bethe-Salpeter equation (1) yields, by contour 
integration in the complex-p plane and the residue theorem, the Salpeter equation (where, 
trivially, p = pi = —p 2 and P = M in the center-of-momentum frame of the bound state) 

Mn\ I M^(Pi)7o[^(p,g)0(g)]7oA2(j>2) 
nP) J (2vr)3 V P - E 1 (p 1 ) - E 2 (p 2 ) 

Ar (gi) 7o [Kip, q) 0(g)] 7o A^(j> 2 ) \ 

P + E^pi) + E 2 (p 2 ) J' { ' 

Multiplying the Salpeter equation (2) from the left by Af(pi) and from the right by Af(p 2 ) 
and using A.f(p) A.f(p) = reveals that any of its solutions 4>(p) will satisfy the constraints 

A+(p0 0(p) A+(p 2 ) = Ar(pi) 0(p) A 2 -(p 2 ) = . 

These constraints halve the number of independent components of any solution 0(p) of the 
Salpeter equation (2); obviously, this has important implications for practical calculations. 
Generally, the Bethe-Salpeter interaction kernel K(p, q) may be decomposed into a sum of 
products of Lorentz-scalar potentials V(p, q) and associated Lorentz structures. Assuming 
the latter to be represented by identical Dirac matrices, generically labelled T, the action of 
the kernel K(p, q) on the Salpeter amplitude 0(g) requested by the Salpeter equation (2) is 



[X(p,q)0(g)] = 5>r(p,<7)r0(q)r 



3 This forms the example par excellence of the trivial po dependence of the propagators requested above. 

4 The assumption of free propagators for the bound-state constituents may, however, encounter a serious 
conceptual problem for the following reason: in quantum field theory the infinite tower of Dyson-Schwinger 
equations relates any n-point Green function to at least one (m ^ n)-point Green function; this entails that 
propagators, being two-point Green functions, and the four-point Green function entering as an interaction 
kernel cannot be viewed as independent. In particular, in quantum chromodynamics (the theory describing 
the strong interactions) free propagators for any of its fundamental quark and gluon degrees of freedom are 
incompatible with the feature of colour confinement, exhibited by this unbroken non-Abelian gauge theory. 
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2.3 Reduced Salpeter Equation 

By the charge-conjugation properties of the energy projection operators Af(p), the second 
term on the right-hand side of the Salpeter equation (2) corresponds to the negative-energy 
components of the Salpeter amplitude 0(p): Aj"(pi) <j)(p) A 2 (p 2 ) = A^(pi) 0(p) [A 2 (p2)] c - 
Assuming that it is justifiable to neglect the contribution of this latter term relative to that 
of the first term on the right-hand side of Eq. (2) leads to the reduced Salpeter equation [5,6] 

[P -E 1 (p 1 )-E 2 (p 2 )} ( f>(p) = J -0^Ai(p 1 ) lo [K(p,q) ( f>(q)} lo A 2 -(p 2 ) , (3) 

forming, unlike Eq. (2), an explicit eigenvalue problem for the bound-state masses P — M. 
Formally, this reduction of the (full) Salpeter equation (2) to the reduced Salpeter equation 
(3) may be effected by imposing as further constraint A^(pi) <f>{p) = or <f>(p) A 2 (p 2 ) = 0; 
as trivial consequence of this, all solutions 4>(p) of the reduced Salpeter equation (3) involve 
only the positive- energy components <f>(p) = A^(pi) <f>(p) A 2 {p 2 ) = A+(pi) <f>(p) [A^(p 2 )] c . 
Physically, this simplification supposes that P —E 1 (p 1 )—E 2 (p 2 ) <C P +E 1 (p 1 )+E 2 (p2) or 
[Po-\-Ei{pi)-\-E 2 {p 2 )\~ 1 <C [P — Ei(pi)— E 2 (p2)]~ 1 holds at the level of expectation values; 
this assumption may be justified for semirelativistic and weakly bound heavy constituents. 
In the center-of- momentum frame of the bound state, P = clearly implies p = p\ = —p 2 . 



2.4 Pseudoscalar Bound States 

To start with, let us focus our interest to those fermion-antifermion bound states which are 
represented by a Salpeter amplitude with the least number of independent components: the 
pseudoscalar states. These are bound states composed of some fermion and its antiparticle, 
thus characterized by mi = m 2 = m and well-defined behaviour under charge conjugation; 
the quantum- number assignment common to all pseudoscalar states is J = for their total 
spin, P = — 1 for their parity, and C — +1 for their charge-conjugation parity: J PC = 0~ + . 
Such states are realized in nature in form of, e.g., the pion as quark-ant iquark bound state. 
Henceforth, all indices i — 1,2 distinguishing the bound-state constituents can be dropped, 
since now, for instance, Ei(p) = E 2 (p) = E(p) = E(p) = ^p 1 + m 2 , withp = |p| = a/p 2 . 

The most general expansion of the Salpeter amplitude 4>(p) over a complete set of Dirac 
matrices would, of course, introduce 16 Salpeter components. However, as a consequence of 
the peculiar projector structure of the Salpeter equation (2), manifesting in the constraints 
mentioned in Sec. 2.2, precisely eight of these components are independent. Specifically, for 
describing pseudoscalar states just two of the latter, labelled <fi(p) and ^(p), are relevant: 



, S H (P) 

v?i(p) p(pT + v?2(p) 



75 



In fact, the above decomposition of <f>{p) applies to all bound states of a spin-| fermion and 
a spin- 1 antifermion for which the quantum number s of the sum of the two fermion spins is 
zero, i.e., to all spin-singlet states carrying s = 0, irrespective of the relative orbital angular 
momentum I of the constituents. Pseudoscalar bound states are just the special case I = 0. 
The one additional constraint truncating the full to the reduced Salpeter equation enforces 
the equality of the Salpeter components (fi(p) and (p 2 (p): flip) = ^(p) = fip)- Thus, all 
spin-singlet, notably all pseudoscalar, solutions of a reduced Salpeter equation (3) will read 

m = <P(P) H{P) J f iP) 75 = 2 fip) A + (p) 75 • 
E[p) 
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2.5 Radial Eigenvalue Equations 

In order to follow the more or less standard route of nonrelativistic reduction, let us assume 
that the Bethe-Salpeter interaction kernel K(p, q) is of convolution type, i.e., is of the form 
K{P, q) — K{p—q), which entails for the potentials Vr(p, q) = Vr(p—q), and that K(p, q) 
respects spherical symmetry, i.e., K(p, q) = K((p—q) 2 ), and thus V T (p, q) = Vr((p—q) 2 ); 
trivially, this means that each momentum-space potential Vr(p, q) is the Fourier transform 
of a spherically symmetric configuration-space potential Vr(r) depending just on the radial 
coordinate r = \x\. In this case, any dependence on the angular variables, encoded either in 
spherical harmonics or vector spherical harmonics, can be split off; this factorization effects 
the reductions of both the Salpeter equation (2) [14,15] and its reduced counterpart (3) [16] 
to equivalent systems of coupled equations for the radial factors of all independent Salpeter 
components. For a particular Lorentz structure of the Bethe-Salpeter kernel (such that the 
index V identifying the Lorentz structure may be suppressed), the interactions experienced 
by bound-state constituents enter the radial eigenvalue equations in form of Fourier-Bessel 
transforms V L (p, q) (L = 0, 1, 2, ... ) of the radial static configuration-space potential V(r), 
defined in terms of spherical Bessel functions of the first kind [17] j n {z), n — 0, ±1, ±2, . . . : 

oo 

V L (p,q) = 8nJdrr 2 j L (pr)j L (qr)V(r), p = \p\ , q = \q\ , L = 0,1,2,.... 
o 

In the following, we intend to exemplify our technique of constructing exact solutions to 
the reduced Salpeter equation by considering only a few of those Lorentz structures of 
the interaction kernel K(p, q) that are commonly used in phenomenological descriptions of 
hadrons, notably, of mesons, as bound states of quarks confined by the strong interactions. 5 



2.5.1 (Full) Salpeter equation 

The Salpeter equation for bound states with spin-sum quantum number s = is equivalent 
to a system of two equations for the two independent radial Salpeter components (p±(p) and 
</j 2 (p)- For pure time-component Lorentz- vector kernels (r<g>r = 7°<g>7°), this system reads 



2E(p) V2 {p) + 



dq q 1 



V (p, q) (p 2 (q) = M(p x {p) , 



2E{p) ( p 1 {p) + 



o 



dg q 2 



Vb(p, q) + —, z —, z Vi(p, q) 



E(p)E(q) 



E(p)E(q) 



yi{q) = M(p 2 (p) . 



5 Since the first encounter [15,18] of unstable solutions of Salpeter equations with confining interactions, 
the stability of solutions of different relativistic bound-state equations has been an issue of concern [19-25] . 
For confining interactions of harmonic-oscillator form in configuration space (in which case all bound-state 
integral equations reduce to easier-to-handlc differential equations in momentum space), the conditions for 
stability of bound-state solutions within the instantaneous Bethe-Salpeter framework (in the sense of their 
spectrum being real, discrete, and bounded from below) have been analyzed for those Lorentz structures of 
the Bethe-Salpeter interaction kernel which are most frequently employed for the description of hadrons as 
bound states of quarks: A rigorous analytic proof of the stability of the bound-state energy spectra entailed 
has been constructed for the reduced Salpeter equation [26-28] and a particular generalization of it [28,29], 
formulated by taking into account the full propagators of the bound-state constituents [12,13,30], obtained 
as the solutions of the Dyson-Schwinger equations for the corresponding two-point Green's functions of the 
bound-state constituents, and this problem has also been discussed for the (full) Salpeter equation [31-33]. 
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2.5.2 Reduced Salpeter equation 

The (compared to the full Salpeter equation) simple and unique energy-projector structure 
of the reduced Salpeter equation (3) guarantees that each spin-singlet bound-state solution 
involves just one independent radial Salpeter component <p(p), since </?i(p) = <f2{p) = <p(p)- 
The system of coupled radial equations thus collapses to a single radial eigenvalue equation 
which reads, for a kernel of time-component Lorentz- vector Dirac structure TcgT = 7°(87°, 



oo 



2E(p) V (p) + - 



1 f dqq 



2 



2 J (2vr) 2 



■2 



Tfl \ p Q 

1+ mw)) VaM+ wmi) VlM . 



f(q) = M ififp) , 



for a kernel of Lorentz- vector Dirac structure TcgT = 7 M <87 M , describing all the interactions 
by an effective vector-boson exchange between the two fermionic bound-state constituents, 



oo 



2 m «p) + J g (2 - j^J) Vo(p, 1) <p{q) — M ,( P ) , 



or, for the linear combination 2 TcgT = 7^*87^+75(875 — 1(8)1 of the kernel's Dirac structure, 

00 

2E(p) V (p) + J ^V (p,q)<p(q) = M V (p) . (4) 


2.6 Special Case: Massless Bound-State Constituents (m = 0) 

In the limit of vanishing masses of the two bound-state constituents, i.e., in the case m — 0, 
the bound-state equations of Sec. 2.5 simplify, of course, still further: the Salpeter equation 
with time-component Lorentz-vector kernel, TcgT = 7°(87°, reduces to the set of equations 

00 

/(lq q^ 
j^y V (p, q) y? 2 (q) = M ip x (p) , 



2 P Mp) + J ^V 1 (p,q)c Pl (q) = Mv 2 (p) , (5) 


whereas the single radial eigenvalue equation emerging from the reduced Salpeter equation 
reads, for interactions with time-component Lorentz-vector Dirac structure TcgD = 7°(8>7°, 

00 

2p V (p)+ 1 - J ^[V (p,q) + V l (p,q)} V (q) = M i f(p) , (6) 


for interactions with Lorentz-vector Dirac structure r <8 T = 7^ (8 7^ , 

00 

2p<p(p) + 2 J ^V ( P , q )<p(q) = M<p(p) , (7) 


or, for kernels of the particularly favourable Dirac structure 2 TcgT = 7 M (87^+75(875 — l<8l, 

00 2 

2p V (p) + J ^V (p,q) V (q) = M V (p) . (8) 


Although we focus to reduced Salpeter equations, we need Eq. (5) at an intermediate stage. 
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3 Configuration- Space Potentials V(r) by Inversion 

All the radial bound-state eigenvalue equations recalled in Sees. 2.5 and 2.6, emerging from 
the reduced Salpeter equation (3) under our reasonable assumption of spherical symmetry, 
are for M ^ 2 E(p) homogeneous linear Fredholm integral equations of the second kind but 
still simple enough that, for sufficiently sophisticated choices of the Salpeter solutions <p(p), 
their underlying configuration-space potentials V(r) can be extracted by analytical means. 
For notational convenience, we define the Fourier-Bessel transforms to configuration space 

oo 

V? L (r) = i L y - J dpp 2 j L (pr) ip(p) , L — 0,1 , 
o 

oo 

T L (r) =i L \llJ dpp 2 j L {pr) E{p) <p(p) , L — 0,1 , 



o 



of both radial Salpeter amplitude ip(p) and free-energy part E(p) ip(p) in momentum space. 
In the following, we illustrate the idea of analytical extraction of configuration-space radial 
potentials V(r) by examining explicit examples: For an appropriate ansatz for the Salpeter 
solution ip(p), by application of the Fourier-Bessel transformation to the momentum-space 
bound-state equation considered we would like to get its configuration- space representation 
in terms of </?l(?") and T L (r); from the latter formulation, we should be able to read off V(r). 
Since in Eqs. (4) and (6)-(8) the potential may absorb any mass M ^ 0, we assume M — 0. 
Of course, the simple procedure sketched above can only be followed if the interaction term 
in the bound-state equation under study contains merely a single Fourier-Bessel transform 
Vl(p, q) of the radial potential V(r), i.e., involves a unique value of L. If, on the other hand, 
both Vq{p, q) and Vi(p, q) enter in the interaction term, as happens, e.g., for any interaction 
of time-component Lorentz- vector structure, a different line of reasoning has to be devised. 
For good reasons, we first invert equations for massless bound-state constituents (Sec. 3.1). 
Then we turn to the more delicate case of bound-state constituents of finite mass (Sec. 3.2). 

3.1 Massless Bound-State Constituents: m = 

In order to be able to deal simultaneously with interaction kernels of Lorentz- vector nature 
and of the (simplifying) linear combination yielding Eq. (4), we introduce a parameter rj by 

2 for T <g> T = 7^ <g> 7^ , 

1 for T <g> T = \ ( 7m <g> 7 M + 7s ® 7s - 1 ® 1) • 
By this definition, the reduced Salpeter equations (7) and (8) may be subsumed in the form 

oo 

/do <7 2 
TT^b Mp, q) <p(q) = m <p(p) . (9) 



7/ 







For this equation of motion, its L = Fourier-Bessel transform is straightforwardly found: 

2 T (r) + r] V(r) y? (r) = M cp (r) . 
Thus, the configuration-space potential associated to mass eigenvalue M = is found to be 

V(r) = - 2 -^f\ . (10) 
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In the case of interaction kernels of the time-component Lorentz- vector Dirac structure, 
i.e., Tcg>r = 7°Cg>7°, the reduced Salpeter equation involves, even for zero-mass bound-state 
constituents, Fourier-Bessel transforms Vl(p, q) of V(r) for more than one L, namely, both 
Vo(p, q) and Vi(p, q). From such type of bound-state equation the potential V(r) cannot be 
recovered by applying to Eq. (6) a Fourier-Bessel transformation of a particular value of L. 
In order to overcome this adverse observation, we recall that, as consequence of the equality 
(fi(p) = v 5 2(p) = vip) °f the two independent radial components of Salpeter amplitudes for 
spin-singlet bound states enforced by the reduced-Salpeter constraint discussed in Sec. 2.3, 
for a definite Lorentz structure TcgF of the interaction kernel the radial eigenvalue equation 
resulting from the reduced Salpeter equation (3) can be found by adding, for ipi(p) = ^(p), 
the two equations that constitute the set of (originally coupled) radial eigenvalue equations 
related to the corresponding Salpeter equation (2); see, for instance, Footnote 2 of Ref. [27]. 
Bearing these findings in mind, we seek, for tentative solutions ip(p) of the reduced Salpeter 
equation (6) for pure time-component Lorentz- vector interaction, the responsible potential 
V(r), with the help of the /w/Z-Salpeter "precursor" (5) of Eq. (6), via a two-step procedure: 

1. After equating the two independent components <pi(p) = ^{p) = fip), we represent 
the decoupled relations arising from the set of equations (5) in configuration space by 
application of the appropriate Fourier-Bessel transformation to each of the relations: 

2To(r) + \/ (r)y9 (r) = M Vo (r) , 

2T 1 (r) + V 1 (r)<p 1 (r)=M<p 1 (r) . 

At this stage we must take into account, by an index L = 0,1, the possibility that the 
potentials Vj,(r) derived from each of these relations by analogy to Eq. (10) can differ: 

2. We assume that the unique configuration-space potential V(r) we are seeking may be 
expressed as a linear combination V(r) = c Vo(r)+Ci V\(r) of the auxiliary functions 
Vo(r) and Vi(r), with yet to be determined (of course, constant) coefficients Co and c\. 
We attempt to find these coefficients by inserting our ansatz for V(r) into the slightly 
more complex and intricate reduced Salpeter equation (6) for M = 0; if we manage to 
deduce thereby a solution for Co and ci, our quest for the potential V(r) is completed. 



3.1.1 Non-Confining Interaction Potentials 

Presumably, the first guess that comes to one's mind for the momentum-space bound-state 
amplitude ip{p) is the exponential, involving a parameter A with dimension of inverse mass: 



(p) = 2A 3/2 exp(-Ap) , A>0, |M| 2 = J dpp 2 \<p(p)\ 2 = 1 . (12) 



For this choice, the L = 0, 1 Fourier-Bessel transforms of <p(p) and m = kinetic term read 

[2 4A 5 / 2 , , [2 4iA 3 / 2 r 

M r ) = v: / 2 , X2N2 » Mn 



2 



n (r 2 + A 2 ) 2 ' riw V 7T (r 2 + A 2 ) 

Tn [2 A ,3/2 3A 2 -r 2 (2 16iA 5 / 2 r 

To(r) = V^ 4A FTW' Tl(r) = V^FTW 
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If the Lorentz structure of the interaction kernel is either a pure vector (T(g>r = 7^®7 M ) or 
a linear combination of vector, pseudoscalar, and scalar (2 r<g>r = 7^(87^+75(8)75 — 1<8>1) as 
summarized by the reduced Salpeter equation (9), our findings (10) immediately entail, for 
a Salpeter amplitude ip(p) of the exponential form (12), in configuration space the potential 

nr) = ^(l-4^yl)> V{0) = -\, V(r) >\. 

77 A \ T z + A z J T] A r->oo Tj A 

In the limit r — > 00, this potential V(r) approaches a finite value. It is thus a representative 
of the class of non-confining interactions, with dependence on the radial variable r (in units 
of A, i.e., for A = 1) for Lorentz- vector Bethe-Salpeter kernels (n = 2) as depicted in Fig. 1. 




Figure 1: Configuration-space non-confining potential V(r) = 1— 4/(r 2 +l), determined by 
inversion of the reduced Salpeter equation (9) with kernels of Lorentz-vector nature (n = 2) 
when assuming an exponential form ip(p) oc exp(— p) of its momentum-space solution <p(p). 
Starting at V(0) = —3, this potential V(r) behaves for r — > 00 like V(r) — > 1 (dashed line). 



For the reduced Salpeter equation (6) with time-component Lorentz-vector Dirac structure 
Tcsr = 7°£g>7° and an exponential ansatz (12) for (p(p), the auxiliary functions (11) become 



X \J r 2 + \ 2 J ) w r 2 + A 2 

Inserting these potentials into Eq. (6) fixes the coefficients c 0j i to Cq — c\ — |, which entails 
a configuration-space potential V(r) with behaviour clearly similar to that shown in Fig. 1: 
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3.1.2 Confining Interaction Potentials 

In order to present also an example for a confining potential, we next consider an amplitude 
of normalized Gaussian form, using a parameter A with dimension of inverse mass squared: 

1/4 °° 

tp{p) = 2^-J exp(-Ap 2 ), A>0, |M| 2 ^|d^ 2 |^)| 2 = l . (13) 

o 

Trivially, the L = Fourier-Bessel transform of this function ip(p) is also of Gaussian form: 

^ o(r)= (^) eXp (~iA)- 

The configuration-space potential entailed by Eq. (9) involves the imaginary error function 
erfi(z) defined in terms of the error function erf (z) [17] by erfi(z) = — i erf(i z); for large r, it 
rises like 16 A 3 / 2 exp(r 2 /4 X)/r) y/n r 4 and thus realizes confinement, as illustrated by Fig. 2: 




1 2 3 4 5 



r 

Figure 2: Configuration-space confining potential V(r) = f| — -) erfi(0 — exp(r 2 /4)/y / 7r, 
found by inversion of the reduced Salpeter equation (9) with Lorentz-vector kernels (j] = 2) 
if trying a Gaussian form <p(p) oc exp(— p 2 ) for its momentum-space solution <p(p). Starting 
at V(0) = —2/y/n ps -1.12837. . . , this potential V(r) behaves like V(r) -> oo for r -> oo. 

A time-component Lorentz-vector Dirac structure can be dealt with along the same lines as 
in Subsec. 3.1.1; however, the resulting expressions are lengthy and not really enlightening. 
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3.2 Bound-State Constituents with Non-Zero Mass: m > 

The extension of the particular inversion technique framed here to the case of nonvanishing 
masses m of the bound-state constituents clearly requires somewhat more careful choices of 
momentum-space amplitudes <p(p). One's first attempt might employ the rational function 

I — 5/2 °° 

*(P) = V£(p^^' ^ >0 ' IMI 2 = /dpP>(p)| 2 = l, (14) 

o 

with a mass parameter /x. Its L = Fourier-Bessel transform is, for /x < m, an exponential: 

<£o(?") = 2 /i 3 / 2 exp(— /x r) , < /x < m . 

For illustrative purposes, we apply this Salpeter amplitude to invert the particularly simple 
reduced Salpeter equation (4) with the linear combination T®T = | (7^(87^+75(875 — 1(8)1) 
of Lorentz structures. In order to evaluate the relevant Fourier-Bessel transform of Eq. (4), 

2T (r) + V(r)(p (r) = M (p (r) V(r) = for M = , (15) 

<A)(r) 

we need, for m/fl, the L = Fourier-Bessel transform T (r) of the kinetic part E(p) (f(p). 



3.2.1 Case < n = m 

For the rational Salpeter amplitude (14), the L = Fourier-Bessel transform of the kinetic 
part is easily pinned down if the parameter /x is chosen to be equal to the constituents' mass 
m; the outcome for T (r) is basically the modified Bessel function K v (z) [17] of order v — 0: 

m , N 8 m 5 / 2 r ^ . , 
T r) = K (mr) . 

71 

Hence, the instant reply of our inversion procedure (15) is the configuration-space potential 
V(r) = — - — K Q (mr) exp(mr) , V(r) > - — ln(mr) , V(r) )• . 

7T r^O 7T r— i-oo 

The monotonic increase of this potential V(r), from its logarithmic singularity at the origin 
r = to its asymptotic value for r — > oo, is depicted for constituent mass m — 1 in Fig. 3. 
The singularity is milder than Coulombic and hence compatible with Herbst's findings [34]. 



3.2.2 Case < fi < m 

Allowing the parameter /x in our ansatz (14) for the Salpeter amplitude <p(p) to be less than 
the constituents' mass m, i.e., if /x < m, performing a contour integration and adopting the 
residue theorem enables us to cast the resulting configuration-space potential into the form 

V(r) = -2 { + ^§= 2 \ - ^ / dppexp[- r (p - „)] J . 

For r — >■ oo, the r-dependent terms of the potential vanish and V(r) approaches a constant: 

\/(r) >■ -2 Vm 2 - /x 2 . 

Examining the r-dependent portion of this potential V(r) for r — >■ by L'Hopital's rule, we 
encounter a logarithmically divergent integral. Hence, since for /x f m the asymptotic value 
of V(r) reduces to zero, the behaviour of V(r) resembles that one found for the case ji — m. 
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Figure 3: Configuration-space non-confining potential V(r) = —(8/tt) K (r) exp(r), found 
by studying the reduced Salpeter equation (4) for bound-state constituents of nonvanishing 
mass m with a kernel of Lorentz structure T&T = | (7^(8)7^+75 £8175— 1®1) when relying on 
a rational function <p(p) oc (p 2 +l)~ 2 as ansatz for the momentum-space solution <p(p). This 
V(r) vanishes in the limit r — > 00 but exhibits a logarithmic singularity at the origin r = 0. 



4 Summary, Conclusions, and Outlook 

In this study, we showed how to establish, by analytical means, exact relationships between 
solutions of Bethe-Salpeter equations and the underlying — instantaneous — interactions: 
for elaborate assumptions about the nature of both interactions and resulting bound states 
culminating in a manageable structure of the equations governing such bound states, this is 
accomplished by postulating particular solutions and reading off the interaction potentials. 
Among others, having these rigorous solutions at one's disposal obviously provides a useful 
or even decisive test when solving Bethe-Salpeter equations numerically by conversion into 
equivalent matrix eigenvalue problems (as, e.g., in Refs. [14-16,35-40]) or when attempting 
to construct approximate models for Bethe-Salpeter solutions, as proposed in Refs. [41,42]. 

Clearly, the three-dimensional Fourier transform of any function which depends only on 
a radial coordinate is just the L = Fourier-Bessel transform of this function. Accordingly, 
due to the simplicity of the bound states inspected some but not all of our reduced Salpeter 
equations, more precisely, Eqs. (4), (7), and (8), are, in fact, equivalent to so-called spinless 
Salpeter equations with, where necessary, appropriately adjusted overall coupling strength 
(for reviews on this latter bound-state equation see, e.g., Refs. [43-45]). Earlier attempts to 
construct exact solutions of spinless Salpeter equations may be found in, e.g., Refs. [46-50]. 
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In order to provide a kind of "proof of feasibility" of the inversion technique constructed 
here, this formalism has been elaborated only for the simplest conceivable problem, that is, 
the one posed by the reduced Salpeter equation. There exist, however, exceptional cases for 
which the above findings apply directly, without changes, also to the full Salpeter equation: 
As recalled in Sec. 2.4, bound states with a vanishing sum of the spins of their constituents, 
such as pseudoscalar states, are represented by only two independent Salpeter components, 
the minimal number of independent components for solutions of the full Salpeter equation. 
Correspondingly, for these states the full Salpeter equation becomes equivalent to a system 
of merely two, in general coupled, equations. For a Lorentz structure of the Bethe-Salpeter 
interaction kernel of, for example, the form 2 T<S>T = 7^ (87^+75 (875 — 1(8)1, one of the latter 
equations does not contain any interactions and is therefore of purely algebraic nature [27]. 
For vanishing bound-state mass, the two equations decouple, and the inversion problem for 
the full Salpeter equation thus becomes identical to that for the reduced Salpeter equation. 

The intention behind this study was to carry out an analysis of purely academic nature. 
Nevertheless, one may ask the legitimate question: to which physical bound states observed 
in nature do the above considerations apply? Section 2.3 confines the validity of the reduced 
Salpeter equation to semirelativistic, weakly bound, heavy constituents; this precludes, for 
instance, the pion but not necessarily pseudoscalar mesons composed of heavy quarks, such 
as i] c and i] b , nor the h c and h b . The range of application to be expected for any full Salpeter 
equation is, of course, much wider. However, a thorough study of the latter equation as well 
as the extraction of a realistic potential are definitely beyond the scope of the present work. 
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